Exceptional Points in Atomic Spectra 
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We report the existence of exceptional points for the hydrogen atom in crossed magnetic and 
electric fields in numerical calculations. The resonances of the system are investigated and it is shown 
how exceptional points can be found by exploiting characteristic properties of the degeneracies, 
which are branch point singularities. A possibility for the observation of exceptional points in an 
experiment with atoms is proposed. 
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The appearance of the coalescence of two eigenstates, 
so-called exceptional points [H, in physical systems de- 
scribed by non-Hermitian matrices has attracted growing 
interest [2, 0, El H @] • Typical systems in which such a 
phenomenon can occur are open quantum systems with 
decaying unbound states. One possibility to describe 
these open quantum systems are non-Hcrmitian Hamil- 
tonians obtained with the complex rotation method Q. 
In this case, the resonances appear as complex eigenval- 
ues whose real and imaginary parts are connected with 
the energy and the resonance width, respectively. The 
eigenvectors have not to be orthogonal in contrast to 
Hermitian Hamiltonians describing bound states in quan- 
tum mechanics. In particular, it is possible that the 
eigenspace for two degenerate eigenvalues is only one- 
dimensional, i.e., there is only one linear independent 
eigenvector. If the system of interest depends on a com- 
plex parameter A (or two real parameters), a branch point 
singularity of two eigenstates can appear at critical pa- 
rameter values A c , which are called exceptional points. 
Exceptional points have been discovered in a broad va- 
riety of physical systems. Among them are acoustical 
systems j^j, atoms in optical lattices 0, Q, and com- 
plex atoms in laser fields [10L Detailed experiments have 
been carried out with resonances in microwave cavities 
@> @, El- However, up to now exceptional points have 
not been found in atoms in static external fields. The 
main reason is that there is only one parameter in the 
cases studied most intensely, viz. atoms either in a mag- 
netic or in an electric field. For the occurrence of ex- 
ceptional points, the parameter space has to be at least 
two-dimensional, i.e., at least two real parameters are 
required, which can be represented by crossed magnetic 
and electric fields. Atoms in static external magnetic 
and electric fields are fundamental physical systems. As 
real quantum systems they are accessible both with ex- 
perimental and theoretical methods and have been used 
for comparisons with semiclassical theories FJ, uM, UA | ■ 



They are ideally suited to study the influence of excep- 
tional points on quantum systems. E.g., the occurrence 
of phenomena like Ericson fluctuations in photoioniza- 
tion spectra has been demonstrated both in numerical 
studies 3, lH and experiments [l7j |. 



In this Letter, we investigate numerically the reso- 
nances of the hydrogen atom in static magnetic and elec- 
tric fields and report the first detection of exceptional 
points in this system. The confirmation of the existence 
of exceptional points supplements the richness of phe- 
nomena which have been found in the spectra of atoms in 
static external fields. Furthermore, we propose a method 
which can be used to verify the existence of exceptional 
points in experiments. 

Exceptional points can occur in systems which are de- 
scribed by a non-Hermitian matrix and depend on one 
complex parameter or two real parameters. At critical 
parameter values A c , the eigenvalues and eigenvectors 
of the matrix pass through branch point singularities 
[H 0, EH- This can easily be understood by studying 
the two-dimensional matrix [l[ 



M(A) = 



1 A 

A -1 



(1) 



with the complex parameter A. The eigenvalues t\ = 
Vl + A 2 and £2 = -yl + A 2 are two branches of one an- 
alytic function in A. An exceptional point occurs, e.g., 
for A = i, where the eigenvectors of the two degenerate 
eigenvalues coalesce and the only linearly independent 
eigenvector reads x(X = i) = (l,i) T . The branch point 
singularity leads to a characteristic behavior of the cor- 
responding eigenvalues under changes of the parameters. 
If one chooses a closed loop in the parameter space and 
calculates the eigenvalues for a set of parameter values 
on this loop, in general, the eigenvalues also traverse a 
closed curve. However, when an exceptional point is en- 
circled, two of the eigenvalues do not traverse a closed 
loop. The two eigenvalues which degenerate at the ex- 
ceptional point are permuted during the traversal of the 
circle in parameter space [1] . This can be seen by plotting 
the paths of the two eigenvalues in the complex energy 
plane. After one circle, the first eigenvalue will travel 
to the starting point of the second and vice versa, as il- 
lustrated in Fig. [T] As a consequence, the path of one 
eigenvalue is not closed, if one traversal of the loop in 
parameter space is performed. But the path is closed if 
the parameter space loop is traversed twice. 

As one would expect, the eigenstates of the system are 
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FIG. 1: (Color online) (a) Circle in parameter space A with 
the exceptional point A c = i as center point for the simple 
model {l}. (b) Eigenvalues £1,2 calculated for the parameter 
values from (a) indicated by squares and diamonds, respec- 
tively. In this special case, each of the two eigenvalues tra- 
verses a semicircle. In (a) and (b) , the filled symbols represent 
the first parameter value Ao and the corresponding eigenval- 
ues ei,2(Ao), respectively. The arrows point in the direction 
of progression. 



permuted in the same way as the corresponding eigenval- 
ues Additionally, after each circle around the excep- 
tional point, one of the two eigenvectors changes its sign 

Q3, e.g., [*i,* 2 ] ^ [*2,-*i]. This phase change 
continues with every additional loop. Therefore, four cir- 
cles are required to return to the initial states. 

In this Letter we investigate exceptional points in the 
hydrogen atom in static external fields. The Hamiltonian 
in atomic units without relativistic corrections and finite 
nuclear mass effects |20j reads 

H= l -p 2 - l -+ l - 1 L z + l - 1 2 (x 2 +y 2 ) + fx 1 (2) 

where L z is the z component of the angular momentum, 
and 7 and / are the field strengths of the magnetic and 
electric fields oriented along the z- and x-axis, respec- 
tively. The only constants of motion are the energy and 
the parity with respect to the (z = 0)-plane. Here, we 
concentrate on states with even z-parity. One possibility 
to calculate the resonances of the system is the complex 
rotation method @, [2l[ . The coordinates of the system 
r are replaced with the complex rotated ones re in 
the Hamiltonian and wave functions. The transforma- 
tion leads to a complex symmetric matrix representation 
of the Hamiltonian. Introducing dilated semiparabolic 
coordinates (see, e.g., [TU) leads to a generalized eigen- 
value problem of the form 



(3) 



with a complex symmetric matrix A("f,f), a real sym- 
metric matrix £?, and the generalized eigenvalue k = 
— (1 + 2b 4 E), where E is the energy and b a complex- 
valued dilation parameter, which includes the complex 
rotation. Above the ionization threshold, resonances are 
uncovered as complex energies E. The computation of 
the eigenvalues was performed by applying the ARPACK 
library [22j to matrices with dimensions of about 10,000 



to 12,300. The algorithm uses the implicitly restarted 
Arnoldi method and solves large scale sparse eigenvalue 
problems efficiently, even for non-Hermitian matrices. In 
general, but not at the exceptional points, the eigen- 
vectors of the resonances can be normalized such that 
(^ilBl^j) = Sij. The external fields 7 and / are the 
two external parameters which determine the eigenstates. 
Exceptional points do exist in atomic spectra if the fields 
can be chosen in such a way that a coalescence of two 
states occurs. The crossed field hydrogen system fulfills 
all necessary conditions for the appearance of exceptional 
points, however, one has to find them in the spectrum to 
prove their existence. 

In practice, it is very difficult to look for exceptional 
points by searching for degeneracies of two complex 
eigenvalues. The variation of the parameter values 7 and 
/ does not lead to clear indications for degeneracies, and, 
it is not known in advance which parameter values arc a 
good choice for starting the search. However, the method 
of encircling a point in the (7, /)-plane and searching for 
eigenvalues which are permuted has been very successful. 
A good choice for such a closed loop is a "circle" in the 
parameter space of the two field strengths with a radius 
a < 1 chosen relative to the center (70, fo): 
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The phase change of one of the eigenfunctions after a 
circle around an exceptional point can be seen as a sign 
change of an arbitrary matrix element pu — (^i | M ^2) ■ 
We chose M = 1, i.e., 



Pij = (*<| 



(5) 



which can be obtained easily in our calculations and is 
not diagonal for the eigenstates which fulfill the orthog- 
onality relation (^i\B\^j) — Sij mentioned above. 

With the method described, exceptional points have 
been found for the first time in the spectrum of the hy- 
drogen atom in static external fields. Fig. [2] (a) shows a 
typical result obtained in a numerical calculation. The 
squares and the diamonds represent one of the eigenval- 
ues at different field strengths, respectively. In this ex- 
ample, using 20 steps on the circle in parameter space 
has been sufficient to obtain a clear signature of the 
branch point singularity. The "radius" of the circle ac- 
cording to Eq. (J4j) was a = 0.01. This value is sufficiently 
large to have an exceptional point inside the circular area 
with a high probability and is small enough to track the 
paths of the eigenvalues with a low number of calcula- 
tions. Fig. [2] (b) shows the position of the degenerate 
eigenvalues (marked with an arrow) in the complex en- 
ergy plane among the resonances in their vicinity. For 
the loop shown in Fig. [5] (a), the matrix element p from 
Eq. (J5J) was calculated. The phase of the complex value, 
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TABLE I: Examples for exceptional points in the spectrum 
of the hydrogen atom in crossed magnetic (7) and electric (/) 
fields. All values in atomic units. 
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degeneracy of the two eigenvalues, the existence of only 
one eigenvector, and the permutation of eigenvalues were 
used to identify the exceptional points. 

In experiments, the complex eigenvalues cannot be ob- 
tained directly. But it is possible to measure the pho- 
toionization cross section and to extract the energy (real 
part of a complex eigenvalue) and width (imaginary part) 
of the resonances. The cross section can be written in the 
form [H 



FIG. 2: (Color online) (a) Paths of the two eigenvalues which 
degenerate at the exceptional point in the complex energy 
plane. The squares and the diamonds represent one of the 
eigenvalues, respectively. Each point of one eigenvalue be- 
longs to a different set of parameters. The path in the field 
strength parameter space is a circle defined in Eq. Q with 
a = 0.01 (see inset). The first set of parameters and the 
corresponding eigenvalues are represented by filled symbols. 
The arrows indicate the direction of progression. The cross 
marks the position of the exceptional point in parameter space 
and the corresponding complex energy of the degenerate res- 
onances, (b) Resonances in the complex energy plane for 
the parameter values 7 = 0.00572 and / = 0.000256 (atomic 
units) at the exceptional point. The position of the two de- 
generate eigenvalues is marked with an arrow. 




FIG. 3: (Color online) Phase of the complex matrix element p 
defined in Eq. <(Sj . It changes its value by ir during the traver- 
sal of the circle in parameter space. The result corresponds 
to the expected sign change. 



plotted in Fig. [3j changes its value by tt. This clearly in- 
dicates a sign change for one of the two eigenstates [l9| 
as mentioned above. 

Examples for parameter values of exceptional points 
and the position of the degenerate eigenvalues in the 
complex energy plane are given in Table U The values 
were obtained by minimizing the distance of two eigen- 
values which indicated a branch point singularity. The 



a(E) = Ana{E - E Q )lm ^ 



*<*>|£>K> 



W\2 



E H - E 



(0) 



where is the initial state with energy Eq, which is 
supposed to be known. The final state with complex en- 
ergy Ej is represented by tyj , D is the dipole operator 
for a given polarization, and a is the fine-structure con- 
stant. The superscript $ on the initial and final states in- 
dicates the angle of the complex rotation used in the cal- 
culation, however, it is important to note that the cross 
section is independent of that angle in converged spec- 
tra. The Fourier transform c{t) (i.e., the time signal) of 
([5]) has the form c(i) = dj ex.p(iEjt) and, therefore, 
the complex energies Ej can be obtained using the har- 
monic inversion method [13] with a high precision. In this 
method, the Fourier transform of the profiles in the pho- 
toionization cross section is used to extract their complex 
amplitudes a,j and energies Ej by solving the nonlinear 
set of equations 



i 



cij expQEjt), 



0,1,. 



,2K 



(7) 



for K values c n = c(t n ) on an equidistant grid. 

With this knowledge, a search for exceptional points 
in an experiment starts with the measurement of the 
photoionization cross section for different field strengths, 
which are located on a closed loop in the (7, /)-space, 
e.g., on a circle of type (j4|). For each measurement, one 
obtains a spectrum which is used to extract the complex 
energies of the resonances. The energy values are drawn 
in a diagram. After these steps, the same method which 
is used to search for branch point singularities in nu- 
merical calculations can be applied to the experimental 
results. The diagrams can be used to look for the char- 
acteristic open curves of single eigenvalues. Fig. Q] (a) 
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FIG. 4: (Color online) (a) Photoionization cross section of 
the energy range in which two eigenvalues connected with an 
exceptional point appear. The cross section is shown for eight 
different pairs of parameter values 7 and / being located on a 
circle around the exceptional point. All values in atomic units, 
(b) Complex energy eigenvalues extracted from the cross sec- 
tions with the harmonic inversion method. Each eigenvalue 
is drawn with a different symbol. Altogether, 16 pairs of pa- 
rameter values were used. In (a), every second cross section is 
shown. The signature of a branch point singularity connected 
with the two eigenvalues labeled with squares and diamonds, 
respectively, is clearly visible. 



shows an example of a typical result for the photoioniza- 
tion cross section for different field strengths which are 
located on a closed loop in the parameter space. The 
loop encircles an exceptional point. Although not ex- 
perimental but theoretical photoionization spectra have 
been used as input for the calculations, it is obvious that 
it is possible to extract the complex energies from a typ- 
ical cross section as it is obtained in an experiment. In 
Fig. [4] (b) one can see the paths of complex eigenvalues 
extracted from the spectra. In our example, 16 cross sec- 
tions were used. This number is sufficient to give a clear 
indication of an exceptional point as shown in Fig. [4] (b). 
The eight spectra in Fig. |4] (a) look rather similar. Note 
that the deviations become more pronounced when the 
parameter space radius is increased. This may be helpful 
for the analysis of experimental spectra with noise. 

In summary, we have found the first exceptional points 
in numerical spectra of an atom in static external fields. 
The branch point singularities can be detected by the 
permutation of two eigenvalues when an exceptional 
point is encircled in parameter space. Further properties 
of the branch point singularities can be used to verify 



their existence. With the harmonic inversion method, 
it is possible to extract the complex eigenvalues of reso- 
nances from experimental photoionization cross sections 
and to detect exceptional points in experimental data. 
This opens the way for the experimental observation of 
exceptional points in atomic spectra. 
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